The Variable Ocean V: Earthquake Precursors

The conventional approach to earthquake warning relies on seismic wave detection. P-wave monitoring systems can provide seconds of alert between the arrival of the faster primary wave and the destructive secondary wave, sufficient time for automated systems to halt trains or open fire station doors, but not a meaningful human warning interval. The question addressed here is whether a genuinely useful precursor signal exists in the hours or days before surface rupture, what its physical basis might be, and whether it is detectable by instrumentation or by biological systems that predate instrumentation by several hundred million years.

Reports of anomalous animal behaviour in the hours before major earthquakes are sufficiently widespread and consistent across cultures and species to merit serious examination rather than dismissal. Dogs exhibiting panic and refusing to remain indoors, birds abandoning roosts simultaneously, snakes emerging from burrows in unusual numbers — these observations span multiple continents and multiple fault systems. A regional study of the Colombian and Venezuelan plate boundary confirmed statistically that domestic dogs are the most consistent species to exhibit pre-seismic behavioural change ahead of major ruptures along that boundary.

The mechanistic question is what these animals are sensing. Three candidate signals present themselves: electrical effects including static charge accumulation and piezoelectric discharge at the surface; infrasound generated by deep crustal strain; and variations in the local geomagnetic field vector. The electrical hypothesis encounters a significant difficulty. Rock has high electrical resistance, and a charge generated at depth would require implausibly large continuous currents to reach the surface intact. More critically, birds observed exhibiting pre-seismic behaviour may be in flight at the time, physically isolated from ground contact and surrounded by air acting as a dielectric. An electrical signal conducted through ground contact cannot account for aerial observations. Infrasound remains a plausible candidate and cannot be dismissed, but it propagates omnidirectionally and does not obviously explain the consistent directional panic response or the apparent ability of animals to orient away from the epicentre. Magnetism survives both objections: magnetic field lines experience no resistance from solid rock and no attenuation from air. A distortion of the local field vector propagates equally to a bird in flight and a dog on the ground.

A study conducted by researchers at the Max Planck Institute provides a more rigorous foundation than anecdotal reports. Rather than relying on human observation, the researchers attached high-precision three-dimensional motion sensors to farm animals — cows, sheep, and dogs — living in a seismically active zone in northern Italy, recording their movements continuously over several months without foreknowledge of when earthquakes would occur. The resulting behavioural dataset was then cross-referenced with a catalogue of thousands of seismic events by independent analysts. The findings warrant attention. The collective behaviour of the animal group showed a statistically significant spike in unusual activity up to twenty hours before earthquakes of magnitude 4.0 and above. The signal was only detectable in the collective: one animal behaving unusually is noise; the entire group shifting its baseline state simultaneously is a measurable pattern. A distance relationship also emerged — animals closer to the future epicentre began reacting earlier, with those directly above it showing anomalous behaviour fifteen to twenty hours in advance, and those fifteen kilometres distant reacting approximately two hours ahead. This linear relationship between proximity and warning time is consistent with a signal that builds from depth and propagates outward through the crust, weakening with distance in a predictable way.

The study additionally noted that the predictive pattern was clearer when animals were confined in barns than when roaming open pasture. The most straightforward explanation for this is not physical amplification by the building structure but a statistical one: confined animals in close proximity to one another produce a cleaner collective signal because their individual movements are constrained and their baseline behaviour is more uniform. Animals spread across open pasture move independently for multiple unrelated reasons, masking the collective response in noise. The barn effect is more likely a measurement condition than a physical phenomenon, and it does not require any additional mechanism beyond the one proposed here.

The Max Planck study has not been independently replicated at scale and its findings should be treated as promising rather than definitive. It is cited here as the strongest available empirical support for the animal precursor observation, subject to the usual proviso that all source claims on this site are treated as unverified until confirmed against primary literature.

In Summary: The animal observation literature points toward a signal that is neither electrical nor acoustic in character. Magnetism is the surviving candidate by elimination rather than by direct proof. The Max Planck farm animal study provides the most methodologically rigorous support to date for a genuine precursor signal, with a measurable collective response and a distance-dependent warning window consistent with a signal propagating from depth. The remainder of this page examines whether a coherent mechanism and instrumental approach can be constructed on that basis.

Magnetic sensitivity is documented across all major vertebrate lineages: mammals, birds, and reptiles. The most parsimonious explanation for its distribution is inheritance from a common ancestor rather than independent evolution in each lineage separately. The most probable candidate is a marine ancestor, for whom three-dimensional magnetic navigation would have had clear survival value. In open water, with no visual landmarks, a compass bearing alone is insufficient; an animal needs azimuth, inclination, and field intensity to establish a three-dimensional positional fix. The cryptochrome proteins implicated in avian magnetic sensitivity — particularly Cry4, expressed in the retina — appear to function by overlaying a magnetic map onto the visual field, processed by a specialised brain region. Human equivalents exist in the form of Cry2 proteins, though conscious access to the resulting signal appears to be suppressed in favour of visual and vestibular navigation.

The implication is that the magnetic sense is not primarily a compass but a three-dimensional positioning system, calibrated to a local field geometry that is normally stable within the timescale of an individual animal's life and movement range. This framing is relevant to the precursor question because it suggests that what triggers the anomalous behaviour is not the detection of an unfamiliar signal so much as the corruption of a familiar reference frame. The animal is not sensing something new; its positional system is receiving contradictory data and generating a threat response. The panic is the output of a navigation system detecting that its map no longer matches reality, not a conscious recognition of seismic danger.

The proposed threshold behaviour follows from this. The animal's magnetic positioning system has a tolerance envelope — a range of field geometries within which it can function and navigate. When the local field vector moves outside that envelope, the system triggers a generalised threat response. The animal does not know why it is panicking. The response scales with the rate and magnitude of field vector change, which means the behavioural signal may intensify as strain builds, giving observers a developing rather than a binary warning.

A large-scale natural experiment supports the proposed link between field vector distortion and animal navigation failure. The South Atlantic Anomaly is a continent-sized region stretching from South America across the South Atlantic toward Africa where the Earth's magnetic field is at its weakest and most geometrically distorted anywhere on the planet. Satellite mapping has established that beneath this region the normal field geometry is severely disrupted, with field lines twisting and plunging toward the core rather than following the smooth inclination angles that characterise the rest of the southern hemisphere. Tracking data for pelagic birds crossing the South Atlantic — Arctic terns, shearwaters, and similar long-distance migrants — shows that flight paths lose their characteristic linear efficiency within this zone, with birds executing wide looping corrections, reversing direction, or making landfall on entirely the wrong continents. The phenomenon is sufficiently well documented to have its own ecological term: vagrancy. Multi-decade tracking datasets have established a statistical relationship between the severity of local geomagnetic distortion and the degree of navigational failure.

This is not an earthquake precursor event. No crustal strain is involved. What the South Atlantic Anomaly demonstrates is simply that when the local field vector departs sufficiently from the geometry an animal's positioning system expects, navigational failure follows. It is the permanent, macro-scale version of the temporary distortion that the proposed earthquake mechanism would produce. The biology works as described: field vector corruption produces navigational breakdown, independently of the cause of that corruption.

The deep origin of the South Atlantic Anomaly is relevant to the broader framework developed on this page. The anomaly sits above the African Large Low-Shear-Velocity Province, a continent-sized body of anomalously dense material at the core-mantle boundary which a number of geophysicists have proposed represents accumulated ancient material that has sunk to the base of the mantle over hundreds of millions of years. The origin of this body is contested in the literature — proposed candidates include ancient subducted oceanic slabs, remnant primordial mantle material, and the reabsorbed remnants of large igneous province events. One possibility worth noting is that the Deccan Traps — the enormous flood basalt province formed approximately 66 million years ago when the Indian subcontinent passed over the Réunion mantle plume — may represent material of this character now resident at the core-mantle boundary beneath the South Atlantic. Geophysicists have traced the root of the Réunion plume to the eastern edge of the African LLSVP, providing a direct mapped connection between the surface volcanic event and the deep mantle structure underlying the anomaly. If a body of that volume and basaltic composition has been progressively absorbed into the mantle and accumulated at the core-mantle boundary, the disruption to core convection dynamics and consequently to the field lines generated there would be proportionally large and permanent. This is offered as a proposed explanatory candidate rather than an established origin, but the Réunion plume track gives it a traceable geological provenance that distinguishes it from the other candidates, and it is consistent with the broader framework developed here in which the Curie release horizon operating on basaltic material at the mantle interface is the fundamental source of magnetic field distortion.

In Summary: Vertebrate magnetic sensitivity is most plausibly a three-dimensional positional inheritance from a marine ancestor. Pre-seismic panic behaviour is proposed to reflect the corruption of a navigational reference frame rather than the detection of danger per se, with panic onset occurring when field vector distortion exceeds the animal's positional tolerance. The South Atlantic Anomaly provides a permanent macro-scale demonstration that field vector corruption produces navigational breakdown in birds, independently of any seismic cause, validating the biological mechanism the earthquake precursor hypothesis depends upon.

Oceanic crust is not uniform in its physical properties. Formed continuously at mid-ocean spreading ridges, it cools progressively as it moves away from the ridge, becoming denser and mechanically stiffer with age and distance. The result is that the lithosphere arriving at a subduction trench is not homogeneous but consists of strips of varying thickness and rigidity reflecting their age at formation.

A further variable is proposed here. The pressure under which oceanic crust forms and cools is not constant but reflects sea level at the time of formation. Higher sea level imposes greater hydrostatic pressure on newly forming crust, producing more effective compression and cooling and consequently denser, more rigid material. Lower sea level produces relatively more ductile crust. Glacial cycles, operating on timescales of roughly 100,000 years for the dominant periodicity, would at typical spreading rates of two to ten centimetres per year produce strips of alternating stiffness on the order of two to ten kilometres in width — a geologically significant scale at the subduction interface. Transitional periods of fluctuating sea level may additionally produce pre-fractured material, weakened by repeated loading and unloading cycles during formation, which would behave differently again from either end member when it eventually reaches the subduction zone: neither maximally stiff nor ductile but internally damaged and liable to fail at lower stress thresholds.

This is not a mechanism that begins with the earthquake and works backward to find a cause. It is a continuous process, analogous in character to the variable ocean itself — a system of differential pressures, densities, and stiffnesses in constant slow interaction, with no discrete start point. The surface rupture that registers as an earthquake is the final and fastest stage of a process that began when the relevant strips first met resistance at the mantle interface.

As this heterogeneous slab descends, adjacent strips of different density and stiffness are forced along the same curved path. Where a stiffer strip meets a more ductile one at depth, differential resistance to bending creates stress concentrations at the boundary between them. Lower in the crust, where temperatures are high enough for plastic deformation, this stress is absorbed gradually. As the strain field propagates upward into cooler and more brittle material, deformation accelerates toward eventual brittle fracture at the surface. The earthquake is not an event; it is a process reaching its terminal stage.

The piezomagnetic effect — the well-documented phenomenon whereby mechanical strain in rock alters the orientation of magnetic mineral domains and thereby distorts the ambient field passing through it — provides the link between this deep process and the surface signal. As strain accumulates at the strip boundaries, the distortion of the local magnetic field vector begins and builds. Because magnetic effects propagate through solid rock without attenuation, this distortion is present at the surface from the earliest stages of the deep process, hours or days before any surface rupture occurs. The signal is not generated by the earthquake; it precedes and predicts it.

A deeper mechanism is proposed to operate at the base of the subducting slab, at the contact boundary between solid basalt and the surrounding mantle, which may account for both the character of the precursor signal and the large-scale anomalies discussed in Section II.

Magnetite and related iron-bearing minerals hold permanent magnetic domain orientations locked in when the rock cooled below the Curie point at the spreading ridge — approximately 580 degrees Celsius for magnetite. Above this threshold, thermal energy overcomes the force holding magnetic domains in alignment and the material becomes magnetically disordered. Mantle temperatures well exceed this threshold. The contact zone between the cold descending slab and the hot surrounding mantle therefore marks a Curie release horizon: a boundary across which previously locked magnetic domains transition from ordered permanent magnetisation to thermal disorder, releasing the field geometry they have been holding into the surrounding medium.

This release is not orderly. The domains being released were locked under strain at angles reflecting the stress history of the rock at the moment of original cooling, subsequently modified by the mechanical history of the slab's journey across the ocean floor. As each layer of the slab crosses the Curie threshold from below, the orientations released into the surrounding field are chaotic and high-amplitude relative to the smooth ambient field. The local field geometry above the contact zone reflects this input continuously.

Critically, the position of the Curie boundary is not fixed. As strain builds in the slab above the contact zone, pressure on the boundary layer changes, altering the local thermal gradient slightly and moving the release horizon up or down through the rock. Each incremental shift releases a new layer of previously locked domains. This is the proposed mechanism for the stepped, non-recovering character of the precursor signal detected at the surface — each increment of the ratchet corresponds to a small migration of the Curie boundary releasing another layer of stressed domains. The signal does not recover between steps because the released domains do not re-lock; the process is thermally irreversible at that depth.

The driving force of the subduction system is not the mantle wedge but the weight of the descending slab itself. As oceanic crust travels away from the spreading ridge it cools, contracts, and becomes progressively denser than the mantle beneath it. At the subduction boundary this dense, cold material begins to sink under its own weight — a process termed slab pull — and once the leading edge has descended sufficiently, gravity continues to draw the remainder of the plate downward. The plate drives the system; the mantle wedge corner flow is a consequence of this descent rather than its cause, generated by the friction of the moving slab dragging adjacent mantle material into rotational circulation. This distinction is important for the precursor mechanism because it means that variations in slab stiffness — the striation thickness variable introduced above — directly affect the friction at the slab-mantle interface and therefore the local thermal gradient at the Curie boundary. A thicker, stiffer striation passing through the contact zone increases frictional resistance, locally elevating temperature at the Curie front and accelerating domain release. The slab is not being processed by the wedge; it is forcing itself through the wedge, and the character of that forcing determines the character of the magnetic signal produced.

The speed at which the mantle wedge corner flow processes arriving striations is not constant but varies with depth and temperature, with direct consequences for the character of seismic output at different subduction zones. At shallower depths the mantle wedge is cooler and more viscous — the rotational flow is sluggish, friction between the descending slab and the wedge material is high, and the rock retains sufficient structural strength to accumulate strain over long periods without releasing it. The Curie boundary migrates slowly through each arriving striation, processing thick strips over extended intervals and generating large, widely spaced magnetic pulses. The seismic output mirrors this: the fault locks over a large area, stress accumulates over centuries, and eventual failure is infrequent and high-magnitude.

At greater depths the ambient mantle temperature is significantly higher and viscosity correspondingly lower. The corner flow vortex rotates faster, shearing action at the interface is more rapid, and the rock at those depths is too ductile to accumulate large strain before releasing it. Striations are processed more quickly regardless of thickness, generating a faster, more continuous sequence of smaller magnetic pulses. The seismic output is frequent, lower-magnitude events rather than periodic catastrophic rupture. Seismic tomography and heat flow measurements across active margins globally place the top of the mechanically coupled wedge — where the descending slab first engages the corner flow — at a consistent depth of approximately 70 to 80 kilometres. Above this depth the plate moves through a relatively stagnant cold forearc region where coupling is minimal. Below it the active rotational processing begins, with its character determined by the local temperature and viscosity of the wedge material at that depth.

Independent evidence for fine-scale magnetic variation within the major seafloor stripes comes from high-resolution deep-tow magnetometer surveys, in which instruments are towed close to the seafloor rather than measured from the ocean surface. These surveys reveal that the large-scale polarity stripes are not internally uniform but contain fine-scale magnetic fluctuations within them, referred to in the literature as tiny wiggles. These variations are attributed to short-period fluctuations in cooling rate and sea level at the spreading ridge — precisely the formation variable proposed in this framework as the source of striation thickness variation. The existence of fine-scale magnetic texture at the ten to twenty metre scale within the major polarity stripes provides direct observational support for the micro-striation sequence the precursor mechanism depends upon, and suggests that the physical barcode the framework proposes is present in the oceanic crust as a measurable feature rather than a theoretical construct.

In Summary: Differential bending between lithospheric strips of varying stiffness, with stiffness partly determined by sea level at time of formation, is proposed as the source of a piezomagnetic precursor signal. At the base of the subducting slab, a Curie release horizon operates continuously as the slab heats through, releasing previously locked magnetic domain orientations into the surrounding field in a stepped, non-recovering sequence. The speed of the mantle wedge corner flow varies with depth and temperature, producing large infrequent events at shallow cool zones and frequent smaller events at deeper hotter zones. The process is continuous from first bending moment at depth to brittle surface fracture, with the magnetic distortion propagating to the surface throughout.

The field distortion predicted by this model is not a simple compass deviation — a shift of north toward east or west — but a tilt of the field vector toward the vertical. Under normal conditions at any given location the field has a characteristic inclination angle, the angle between the field vector and the horizontal plane, which varies with latitude. The proposed precursor signal is a slow drift of this inclination away from its established baseline, either as a steady progression or as an accumulation of small incremental shifts that fail to fully recover between episodes, producing a ratcheting pattern that reflects the stepped nature of the bending process as successive strip boundaries reach their stress threshold.

This signature has a useful property: it differs in character from the principal sources of geomagnetic noise. Solar wind variation and ionospheric disturbance predominantly affect the horizontal components of the field and produce oscillatory signals that return to baseline. A tectonic inclination drift would be directional and non-recovering within the precursor window. The vertical component of the field is additionally the quietest channel in terms of external electromagnetic interference, making it potentially the most favourable for extracting a weak tectonic signal against background noise.

The precursor signal measured at the surface is proposed to be not merely an indicator of strain accumulation but a direct playback of the magnetic record encoded in the seafloor striations as they are sequentially erased at the Curie boundary. When oceanic crust formed at the spreading ridge, each striation locked in the orientation of the ambient magnetic field at the time of cooling — alternating between normal and reversed polarity across successive strips, with amplitude varying according to the iron content and volume of material in each strip. This record has been preserved in the rock for tens of millions of years. As each striation reaches the Curie boundary at depth and crosses the thermal threshold, its locked domain orientations are released and its contribution to the local field is abruptly removed. The surface sensor does not detect new energy being generated; it detects the sequential erasure of an ancient record, striation by striation, as the plate passes through the thermal front. The frequency of the resulting surface signal reflects how rapidly successive striations are being processed — a function of plate speed and striation width. The amplitude reflects the volume and iron content of each striation — a function of its formation history and sea level at time of deposition. The precursor signal is in this sense a real-time broadcast of the plate's magnetic biography, transmitted to the surface at the moment of its destruction.

Pre-seismic ultra-low frequency electromagnetic anomalies recorded before major subduction zone earthquakes provide instrumental corroboration for the playback mechanism proposed here. In the period before the 2011 Tohoku earthquake in Japan, ground-based magnetometers recorded slow pulsating shifts in the local magnetic field vector in the days preceding rupture. Under conventional seismic models these signals are difficult to account for because no mechanical fracture had yet occurred at the surface. Under the framework proposed here they are the expected surface expression of thermal demagnetisation at the Curie boundary — the sequential erasure of micro-striation packets at depth producing a slow, pulsating field disturbance that propagates to the surface before the locked zone above reaches mechanical failure. The ULF signal is not a mystery requiring a separate explanation; it is the precursor signal this framework predicts, observed instrumentally before one of the largest recorded subduction zone earthquakes.

In Summary: The predicted precursor signal is an inclination drift toward vertical — slow, directional, and non-recovering — which differs in statistical character from solar and ionospheric noise and occupies the quietest portion of the geomagnetic spectrum.

Existing geomagnetic monitoring networks, most notably INTERMAGNET, record field variations continuously using three-axis fluxgate magnetometers. The precursor signal has not been reliably extracted from this data, though a number of studies have reported field anomalies in the 24 to 72 hour window before major events. The difficulty is signal-to-noise ratio: measuring the full field and looking for small deviations against a background of solar and ionospheric variation is a demanding extraction problem even with modern signal processing.

A geometrically motivated alternative is proposed. If a fluxgate sensor is oriented not along the ambient field baseline but perpendicular to it — specifically, oriented to measure in the vertical plane while the normal field at that location is predominantly horizontal — the sensor sits in a mechanical null with respect to the background field. The normal diurnal variation, predominantly horizontal in character, exerts minimal influence on a sensor oriented perpendicular to it. The predicted precursor signal, a tilt of the field toward vertical, would by contrast strike such a sensor at an increasingly direct angle as strain builds, producing a clean positive deflection from null rather than a small deviation from a large background value. This geometry pre-subtracts the dominant noise source in hardware rather than software, which may offer advantages in stability and drift characteristics over purely computational baseline subtraction.

Slow electrical interference from motor noise and similar sources, which alternates at 50 to 60 Hz, is readily removed by a low-pass filter either in hardware or software, leaving only the slow multi-hour tectonic drift exposed. The sensor array does not require artificial intelligence or complex pattern recognition to extract the signal; the geometry of the instrument does the filtering work before the data reaches a computer.

The instrumental prediction of this framework is specific: a sensor array of this type deployed near active subduction zones should record slow inclination drift, potentially with a ratcheting character, in the hours to days before significant rupture events. This is testable against existing earthquake catalogues if the sensor geometry can be replicated in analysis of historical INTERMAGNET data from stations near well-documented events.

In Summary: A fluxgate sensor oriented perpendicular to the local field baseline would sit at mechanical null against normal diurnal variation while registering the predicted inclination drift as a clean positive signal. The geometry is the filter. This is a testable instrumental proposition.

If the proposed mechanism is valid, the precursor signal may carry information beyond the simple fact that strain is accumulating. This is not presented as a certainty but as a logical consequence of the framework developed above, offered as a set of testable propositions.

The plate speed and earthquake recurrence interval for Cascadia provide a concrete numerical check on the framework's consistency. The Juan de Fuca plate subducts at approximately four centimetres per year. Over a 500-year recurrence interval — the documented average for major Cascadia ruptures — this produces 20 metres of accumulated plate displacement. Paleotsunami modelling of the last major Cascadia rupture in January 1700 independently estimates a fault slip of 15 to 20 metres, consistent with that accumulated displacement to within measurement uncertainty. The framework predicts that this 20-metre cycle corresponds to the processing of a specific packet of micro-striations through the Curie boundary — the thermal playback of approximately 20 metres of seafloor magnetic record driving the friction and stress accumulation that culminates in rupture. The agreement between the plate speed calculation, the recurrence interval, and the independently estimated fault slip is consistent with the proposed mechanism, though it does not uniquely confirm it. It is noted here as a numerical check rather than proof, and would require replication across multiple subduction zones with differing plate speeds and recurrence intervals before it could be treated as strong evidential support.

A coastal sensor array measuring inclination drift across multiple stations would potentially yield three independent variables, each suggesting something different about the developing event beneath.

The spatial extent of the anomaly across the array would correspond approximately to the along-strike length of the stressed slab section. Rupture length is a primary determinant of earthquake magnitude — the 2004 Sumatra event involved approximately 1,200 kilometres of rupture along the Sunda subduction zone, which is why it produced a magnitude 9.1 rather than a 7. If sensors across a 400-kilometre coastal span all begin recording inclination drift within the same time window, that spatial coherence suggests the stressed section is of comparable scale. A signal confined to a 50-kilometre span suggests something considerably smaller. The array does not predict magnitude directly, but it maps the geometry of the stress zone, from which magnitude range follows as a reasonable inference.

The angle of attack of the subducting plate, developed as a framework elsewhere on this site, interacts with the depth-viscosity relationship to determine the maximum possible magnitude at any given subduction zone, and real-world seismicity confirms the predicted pattern. At shallow subduction angles of ten to thirty degrees the descending plate scrapes horizontally beneath the overriding plate for considerable distance before reaching the 70 to 80 kilometre coupling depth. The contact zone between the two plates within the cold, high-viscosity upper mantle is therefore extensive, friction is maximised across a large area, and the fault can remain locked for centuries accumulating strain. Cascadia and southern Chile exemplify this configuration — shallow angle, massive overriding continental plate acting as a rigid buttress, extended cold contact zone, and a seismic record characterised by infrequent megathrust events in the magnitude 9 range with long intervening periods of near-total seismic silence. The sensor array prediction for such a zone would be large, widely spaced precursor pulses building over an extended period before a single high-magnitude release.

At steep subduction angles of sixty degrees or more the descending plate passes quickly through the cold upper zone and enters hotter, lower-viscosity mantle almost immediately. The cold contact zone is minimal, the fault cannot accumulate large strain before releasing it, and the seismic output is frequent lower-magnitude events. The Mariana Trench exemplifies this configuration — near-vertical subduction, rapid entry into hot ductile mantle, and no recorded megathrust event in the historical record despite continuous seismic activity. The sensor array prediction for such a zone would be a rapid, lower-amplitude precursor hum rather than large stepped pulses. The coupling depth measurement of 70 to 80 kilometres provides an empirical anchor for this distinction — it is the point at which plate geometry and local mantle temperature combine to determine which of these regimes will dominate. The geometry of the array response across these different configurations may therefore suggest something about the plate configuration driving the event, independent of other measurements.

The depth at which the field is being distorted offers a third inference. A magnetic field bent at great depth produces a more diffuse surface signature than one bent in the shallow crust, because the distortion has propagated further before reaching the sensor. If the inclination drift is small in absolute terms but spatially coherent across a wide array, the implication is that the source is deep and the energy budget is large — a small signal from far down represents considerably more stored energy than the same signal originating shallow. Conversely, a large and rapidly developing drift confined to a small area suggests a shallower, more localised event approaching failure quickly. The combination of drift magnitude, spatial extent, and rate of development may together allow a rough depth estimate and consequently a severity range.

Taken together, these three variables — spatial array extent, drift signature character reflecting angle of attack, and inferred source depth from drift magnitude and coherence — suggest that a well-deployed coastal sensor network would produce not a binary warning but a developing picture: an early probability range for event magnitude and affected zone, narrowing as the signal evolves, in the same way that a meteorological forecast narrows from a probability distribution toward a specific prediction as more data arrives. This is offered as a logical consequence of the framework rather than a proven capability, and would require systematic testing against historical earthquake sequences before any operational confidence could be justified.

The angle of attack and plate speed variables together produce a two-axis predictive framework for subduction zone behaviour that is worth stating explicitly as a testable consequence of the mechanism. The angle of attack determines how much of the descending plate remains in contact with the cold, brittle upper mantle before reaching the coupling depth — and therefore how large a contact zone is available to accumulate strain. A shallow angle maximises this contact zone and consequently the potential magnitude of eventual rupture. Plate speed determines how rapidly successive striations are fed through the Curie boundary — and therefore the frequency at which strain is loaded and released. A fast-moving plate processes striations rapidly, resets the stress cycle quickly, and produces frequent events. A slow-moving plate accumulates strain over centuries before a striation of sufficient stiffness finally tips the system into failure.

Combining these two variables produces a predictive matrix consistent with observed global seismicity. A shallow angle with slow plate speed — the configuration of Cascadia — produces the highest possible magnitude at the lowest frequency: strain accumulates over centuries across a vast cold contact zone before infrequent catastrophic release. A shallow angle with faster plate speed — southern Chile — produces similarly high magnitudes but more frequently, as the thermal playback rate is higher and the system resets more rapidly. A steep angle with fast plate speed — the Mariana configuration — produces frequent low-magnitude events as the plate bypasses the cold contact zone almost immediately and enters ductile mantle where strain cannot accumulate to catastrophic levels. This matrix does not require new measurement infrastructure to apply; plate speed and subduction angle are already measured parameters at all active margins globally, and their combination under this framework produces magnitude and frequency predictions testable against the historical record.

The striation thickness variable introduced in Section III has a direct implication for the character of the precursor signal and the irregularity of the earthquake record. Each striation arriving at the subduction interface has a different thickness reflecting its sea level formation history. Thickness determines how long the slab takes to heat through to the Curie release threshold. A thick strip — formed under high sea level pressure, dense and stiff — presents a larger volume of locked magnetic domains to the contact zone and takes longer for the Curie boundary to migrate through it. When that boundary reaches a critical depth within the strip, the volume of domains crossing the threshold simultaneously is large and the resulting pulse in the surface field is correspondingly strong. A thin strip heats through faster, releases its domains over a shorter interval, and produces a smaller, faster pulse.

The sequence of striations arriving at any given subduction zone is determined by the sea level history encoded in the ocean floor at the spreading ridge, which was irregular — glacial cycles of varying intensity, duration, and transition speed producing strips of varying thickness in an irregular sequence. The result at the subduction interface is an irregular series of domain releases: varying pulse magnitude, varying pulse duration, varying interval between pulses. This irregularity in the deep magnetic pulse sequence is proposed as a contributing factor to the observed irregularity of earthquake timing and magnitude at subduction zones. No two subduction zones consume identical sea level histories, which would account for why each zone has a characteristic but non-repeating seismic signature.

For the surface sensor array this has a practical implication. The precursor signal will not be a smooth inclination drift but a stepped sequence of pulses, each corresponding to a striation boundary crossing the Curie release threshold. The rate of stepping and the amplitude of each step may carry information about the thickness and stiffness of the material currently at the interface and consequently about the scale of the impending surface event. A sequence of large, widely spaced steps suggests thick, stiff material building toward a high-magnitude release. Rapid small steps suggest thinner, more ductile material releasing energy incrementally rather than accumulating it toward a single large event. This is offered as a logical consequence of the framework rather than an established result.

In Summary: Spatial array extent suggests rupture length and therefore magnitude range. Drift signature character reflects angle of attack geometry. Drift magnitude relative to spatial coherence implies source depth and energy budget. The stepped character of the precursor signal reflects the striation thickness sequence arriving at the subduction interface, with large widely-spaced steps implying high-magnitude accumulation and rapid small steps implying incremental release. Together these variables suggest a developing forecast rather than a binary prediction, testable against historical data before any operational deployment.

The mechanism proposed in this page does not stand in isolation from the broader body of work on this site. The Variable Ocean model and the Cumulative Thermal Lag framework developed elsewhere treat sea level history as a measurable record legible in the landscape, in place name distributions, and in the coastal and tidal archaeology of southern Britain. The proposal here that sea level at time of crust formation is an input variable into lithospheric stiffness represents a potential extension of that framework into deep time and into the mechanics of the ocean floor itself.

The connection is speculative and is not presented as established. What it suggests is that the same variable — sea level — which shaped the human landscape of the last two millennia may also have encoded structural properties into oceanic crust formed over the last tens of millions of years, properties that influence where and how severely the Earth's surface fractures today. If that connection proves to have merit, sea level reconstruction becomes relevant not only to coastal archaeology and landscape interpretation but to long-term seismic hazard assessment at subduction zones. That is a testable proposition, and it is offered here in that spirit.

In Summary: The sea level variable that anchors the Variable Ocean and CTL frameworks on this site is proposed as an input into lithospheric stiffness at formation, potentially linking coastal landscape archaeology to deep-time subduction mechanics. The connection is speculative and offered for testing rather than as an established result.

The Curie release mechanism proposed at active subduction zones may operate at a vastly larger scale where ancient basaltic material has accumulated at the core-mantle boundary over geological time. This section is the most speculative in the page and is flagged explicitly as such. It is offered as a logical extension of the proposed mechanism rather than as a claim with direct evidential support.

The active mechanism at subduction zones is the mantle wedge corner flow — a rotational circulation of mantle material driven by the descending slab, which brings intense heat directly to the contact zone and acts as the grinding interface through which arriving striations pass the Curie threshold. It is this rotational flow pressing hot mantle material against the cold incoming slab that drives the domain release described in Section III, processing successive striations as they arrive and generating the stepped precursor signal detectable at the surface. The mantle wedge temperature at this interface is in the region of 1300 degrees Celsius — sufficient to cross the Curie threshold, but moderate enough that striation thickness remains the rate-limiting factor, with thicker strips taking longer to heat through and releasing their domains in larger, more widely spaced pulses.

The African Large Low-Shear-Velocity Province is a continent-sized body of anomalously dense material mapped at the core-mantle boundary beneath Africa and the South Atlantic. Its origin is contested. One candidate, noted here because it carries a traceable geological provenance, is the Deccan Traps flood basalt — an enormous volume of basaltic material extruded approximately 66 million years ago when the Indian subcontinent passed over the Réunion mantle plume. Geophysicists have traced the root of the Réunion plume to the eastern edge of the African LLSVP, providing a direct mapped connection between the surface volcanic event and the deep mantle structure underlying the anomaly.

At this depth and scale the same rotational flow principle applies, but the grinding interface is the convecting outer core pressing against the base of the blob rather than a mantle wedge acting on a descending slab. The temperature at the core-mantle boundary is in the region of 3500 to 4000 degrees Celsius — far exceeding the subduction zone mantle wedge and imposing a far steeper thermal gradient on the basaltic material at the contact surface. At this temperature the Curie threshold is crossed rapidly regardless of material thickness, and the rate-limiting factor shifts from striation thickness to the rate at which the outer core's rotational flow can press fresh material against the boundary. Domain release at the contact surface is therefore faster and more continuous than at the subduction interface, producing a steady high-amplitude field disturbance rather than the stepped pulse sequence the subduction mechanism generates. This explains why the SAA presents at the surface as a broad stable anomaly rather than a pulsating signal — the extreme temperature overwhelms the striation thickness variable that produces irregularity at shallower depths.

Rock conducts heat slowly, and the interior of the blob has not yet reached the Curie threshold — locked domain orientations from the original formation and subduction history of the material remain intact at the centre. The Curie boundary is therefore not a surface but a shell, migrating inward through the blob over millions of years as heat penetrates from the outside. Only the outermost layer in direct contact with the core heat source has crossed the threshold and is actively releasing its locked domains. The interior remains magnetically structured and cold relative to the surrounding environment.

The blob cannot descend further into the liquid outer core because liquid iron at that depth is significantly denser than the silicate rock of the blob. The descending material therefore flattens against the core-mantle boundary, spreading laterally and maximising its contact surface with the liquid iron beneath. This geometry is significant: rather than a point contact, the blob presents a broad flat base to the thermal source, ensuring that the Curie release horizon operates across the full lateral extent of the accumulated material simultaneously. The relatively sharp edges of the South Atlantic Anomaly at the surface correspond to the edges of this flattened contact zone — the anomaly boundary marks where the blob's base meets the core, and beyond those edges the normal field geometry resumes. As the inward-migrating Curie shell consumes more of the blob's interior over geological time, the surface expression of the anomaly will evolve, its character changing as the ratio of processed to unprocessed material shifts. It is not a permanent feature of the planet but a stage in the thermal processing of a finite ancient basaltic body, large-scale and non-recovering until the whole intrusion has been processed through the thermal threshold and its locked domain orientations fully released.

The scaling argument developed in this section suggests a further consequence at the extreme end of the process. As the basaltic body at the core-mantle boundary grows through continued subduction input over geological time, the thermal insulation it presents to the liquid iron beneath it increases proportionally. The disruption to core convection — already sufficient to produce the reverse flux patches observed beneath the African LLSVP — would intensify as the contact area expands and the inward-migrating Curie shell processes more of the interior. The reverse flux patches documented beneath the SAA are localised zones where the magnetic field lines point in the opposite direction to the primary dipole. If the volume of material producing this effect grows to the point where the reverse flux energy beneath the blob exceeds the organising force of the primary convection columns, the global dipole loses coherence. The field passes through a period of weakened, multi-polar character before reorganising in the orientation dictated by whichever convection pattern re-establishes dominance. The turbulent circulation forced around the edges of the blob by this obstruction generates, by magnetohydrodynamic induction, reversed field components in the liquid iron — the disruption does not merely weaken the primary dipole but actively produces a competing reversed field whose strength scales with the intensity of the turbulent flow. This is the mechanism of a geomagnetic reversal as observed in the palaeomagnetic record — not a sudden external event but the terminal consequence of sufficient basaltic accumulation at the core-mantle boundary disrupting the dynamo from below. The cause is the same process operating at this page's smallest scale; only the volume of material and the depth of processing differ.